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Abstract 

We present the mathematical background of a software package 
that computes triangulations of mapping tori of surface homeomor- 
phisms, suitable for Jeff Weeks's program SnapPea. The package is 



an extension of the software described in [BriOO]. It consists of two 
programs, jmt computes triangulations and prints them in a human- 
readable format, jsnap converts this format into SnapPea's triangu- 
lation file format and may be of independent interest because it allows 
for quick and easy generation of input for SnapPea. As an application, 
we obtain a new solution to the restricted conjugacy problem in the 
mapping class group. 



1 Introduction 



In [ BriOO ], the first author described a software package that provides an 
environment for computer experiments with automorphisms of surfaces with 
one puncture. The purpose of this paper is to present the mathematical 
background of an extension of this package that computes triangulations of 
mapping tori of such homeomorphisms, suitable for further analysis with Jeff 
Weeks's program SnapPea |[MW92|| .F1 

*This research was partially conducted by the first author for the Clay Mathematics 
Institute. 
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Pseudo-Anosov homeomorphisms are of particular interest because their 
mapping tori are hyperbolic 3-manifolds of finite volume |[Thu98|| . The soft- 



ware described in |BriOO|| recognizes pseudo-Anosov homeomorphisms. Com- 
bining this with the programs discussed here, we obtain a powerful tool for 
generating and analyzing large numbers of hyperbolic 3-manifolds. 

The software package described in [ BriO0|| takes an automorphism <fi of 
a surface S with one puncture (given as a sequence of Dehn twists) and 
computes the induced outer automorphism of the fundamental group of S, 
represented by a homotopy equivalence / : G — > G of a finite graph G C S 
homotopy equivalent to S, together with a loop a in G homotopic to a loop 
around the puncture of S. The map / and the loop a determine up to 
isotopy |BriO(J| , Section 5.1]. 

In Section |2], we describe an effective algorithm for computing a trian- 
gulation of the mapping torus of 0: S — > S, given only / : G — » G and a 
(Theorem |2.3j ). We also present an analysis of the complexity of this algo- 
rithm (Proposition |2.4|) . The first part of the software package is a program 
(called jmt) that implements this procedure. The program jmt prints its 
output in an intermediate human-readable format. 

In Section ^ we explain how to use the software discussed here and the 
isometry checker of SnapPea to solve the restricted conjugacy problem in 
the mapping class group (i.e., the question of whether two pseudo-Anosov 
homeomorphisms are conjugate in the mapping class group). This problem 
was previously solved in [|Mos86|| and ||Hem79|| . One distinguishing feature of 



our solution is that much of it has already been implemented. 

Appendix [A] discusses the second program in the software package (called 
jsnap), which converts the intermediate format of jmt into SnapPea's trian- 
gulation file format. Since SnapPea's format is rather complicated, it is not 
easy to generate input files for SnapPea, and jsnap may be of independent 
interest because it allows users to generate input for SnapPea without having 
to understand SnapPea's file format. 

Finally, in Appendix ||, we present some sample computations that ex- 
hibit some of the capabilities of the combination of SnapPea and the software 
discussed here. 

Immediate applications of the software described here include an exper- 
imental investigation of possible relationships between dynamical properties 
of pseudo-Anosov homeomorphisms (as computed by the first author's train 
track software) and topological properties of their mapping tori (as computed 
by SnapPea). For example, one might look for a relationship between growth 
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rate and volume. Another area where the package described in this paper 
has already been used is the study of slalom knots as introduced by Norbert 
A'Campo ||A'C98|| . 

The software package is written in Java and should be universally portable. 
The programs jmt and jsnap are command line software and can be used 
to examine a large number of examples as a batch job. A graphical user 
interface with an online help feature is also available. 

The package, including binary files, source code, complete online docu- 
mentation, and a user manual, is available at |http : / /www . math . uiuc . edu/ 
"brinkman/. 

We would like to thank Mladen Bestvina and John Stallings for many 
helpful discussions, as well as Jeff Weeks and Bill Floyd for explaining Snap- 
Pea's intricacies. We would also like to express our gratitude to Kai-Uwe 
Bux for critiquing an early version of this paper. 

2 Computing triangulations 

Let (ft : S — ► S be an automorphism of a surface S with one puncture, repre- 
sented by a homotopy equivalence / : G — > G of a finite graph G and a loop 
a in G representing a loop around the puncture of S (see Section [I]). There 
is no loss in assuming that / : G — > G maps vertices to vertices and that the 
restriction of / to the interior of each edge of G is an immersion. 

In this section, we outline an effective procedure that computes a trian- 
gulation of the mapping torus of given only / and a. To this end, we 
construct a simplicial 2-complex K and a face pairing e with the following 
properties. 

1. The space \K\ is homeomorphic to a torus. 

2. For each 2-simplex A of K, there exists a 2-simplex A' of K and an 
orientation reversing simplicial homeomorphism e A : A — > A' such that 

e A' = e A- 

3. The space K/e is homotopy equivalent to the mapping torus of /. 

4. If we let M = (cone over K)/e and obtain M' from M by removing the 
cone point, then M' is a 3- manifold (in particular, the links of vertices 
in M' are 2-spheres). 
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In this situation, M' is homotopy equivalent to the mapping torus of /, 
which in turn is homotopy equivalent to the mapping torus of 0. As M' 
is a 3-manifold, M' is homeomorphic to ||Joh79| , page 6]. 

The triangulation of K induces a triangulation of M, i.e., the tetrahedra 
of M are cones over the triangles of K. The vertices of K give rise to finite 
vertices of M, and the cone point is an ideal vertex corresponding to the 
torus cusp of M^. By computing the links of vertices, SnapPea recognizes 
finite vertices (whose links are 2-spheres) and ideal vertices (whose links are 
tori or Klein bottles). 

Hence, we have reduced to problem of constructing a triangulation of the 
mapping torus of to the construction of the 2-complex K and face pairing 
e, given only the homotopy equivalence / : G — > G and the loop a. The 
construction of K and e is the purpose of the remainder of this section. 

The construction of K and e proceeds in two steps. We construct the 2- 
torus T by gluing annuli using Stallings's folding construction ||Sta83|| . Then 



we construct a triangulation and a face pairing for each of the annuli. 



2.1 Step 1: Subdividing and folding 



We review the notion of subdividing and folding [(3ta83| , |BH92|| . Let G, G' 



be finite graphs, and let / : G' — > G be a map that maps vertices to vertices 
and edges to edge paths. 

If / fails to be an immersion, then there exist two distinct edges a, b in 
G' emanating from the same vertex such that f(a) and f(b) have a nontrivial 
initial path in common. We construct a new graph G[ by subdividing a (resp. 
b) into two edges a±, a 2 (resp. bi, b 2 ). 

Now / factors through G[, i.e., there are maps s : G' — » G' x and g : G' x — > 
G such that / = gos. Moreover, we can choose s and g such that s(a) = aia 2 , 
s(b) = &1&2 and g(ai) = g(bi). We obtain a new graph G' 2 from G[ by 
identifying the edges at and b\. Then g factors through G' 2 , i.e., there is 
a map h: G 2 — > G such that g = h o p, where p is the natural projection 
p: G[ — > G 2 (see Figure |T]). We refer to this process as folding a x and b\. 

Remark 2.1. The notion of folds used in pH92|| differs slightly from that 



introduced in [5ta83 . In [BH92j| , the authors consider homotopy equivalences 
/ : G —>■ G, and folding changes both the domain and the range of /, whereas 
in |Sta83|1 , the author considers maps / : G' — ► G, and folding only affects 



the domain G' . The notion of folds used in this paper is a slight modification 
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Figure 1: Subdividing and folding. 



of the folds in |[5ta83 



The remainder of this subsection details how to construct a sequence of 
graphs and maps 

s~1 y^i SO ^ VO y^l S n _l ™ Pn — 1 ^ 9n 

U — Lto > Lti > Lr 2 • • • L*2n-2 ¥ ( - r 2n-l ¥ ( - T 2n ¥ 

such that 

f = 9n° Pn-1 ° S n _l O . . . O p O S , 

where Si : G 2i — > G 2i+ i is a subdivision, p^ : G 2i+ i - ► G 2i+2 is a Stallings fold, 
and g n : G 2n — ► Co is an immersion. Since / is a homotopy equivalence, g n 
will be onto, hence a homeomorphism. Moreover, for each i = 0, . . . , 2n, we 
will construct a loop ex, in Gj corresponding to a loop around the puncture 
of S. 

Let /: G — > G be induced by a homeomorphism 0: S S, and let a 
denote an edge loop in G corresponding to a loop around the puncture of S. 
Let G = Go, go = f : Go — > G, and ct = a. 

Suppose that go is not an immersion. Then there exist two edges a, b 
emanating from the same vertex in Go such that go(a) and go{b) have a 
common initial segment. Since go is induced by the homeomorphism <ft: S — > 
S 1 , we can find a and b such that a and b are adjacent in the embedding of 
Go in S. 

Since the loop do in G is homotopic to a loop around the puncture, a 
and b will be adjacent in the spelling of ctq. Hence, we can detect a and b 
algorithmically by looking for cancellation between the images of adjacent 
edges in the spelling of oq. 

We obtain G\ from Go by subdividing a and b, and we obtain G 2 from 
Gi by folding the initial segments of a and b. As above, we construct maps 
s : G — > Gi, po~- Gi — > G 2 , and G 2 — > G such that g — Pi ° Po ° s o- 
Let o"i = s (cr ) and obtain cr 2 from po( cr i) by tightening. Since the edges a 
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Figure 2: The graph G embedded in the surface S. The corners of the 
octagon correspond to the puncture of S. Two faces of the octagon are glued 
via an orientation-reversing map if the edge labels match up. 

and b are adjacent in the embedding of G in S, the embedding of Go m S 
induces an embedding of G\ and Gi in S, and o\ and a% are homotopic to 
cr in S. 

Note that the size of g%, i.e., the sum of the lengths of the images under 
gi of the edges in G2, is strictly smaller than the size of go- Hence, after 
repeating this construction finitely many times, we reach a map g n : G2 n —* G 
that cannot be folded and thus has to be an immersion. We have found the 
desired sequence of subdivisions and folds. 

Example 2.2. Let G be the graph with one vertex and four edges, labeled 
a, . . . , d, embedded in a punctured surface S of genus 2 as shown in Figure 
We consider the following homotopy equivalence / : G —> G induced by an 
automorphism of S. 



f(b) 

m 



acdcb 



bcdcbcdcb 



bcdd 



ddcbd 



a 



ababcdcd 
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In f(a), cancellation occurs between the underlined parts of f(a) and /(&), 
and we subdivide a and b in preparation for folding, which gives us the maps 
so: G — > Gi and g\\ G\ — > G (see Figure [[]). In order to reduce notational 
complexity, we only change the labels of those edges that are subdivided. 



so (a) 


= a\a2 


s (6) 


= hb 2 


so(c) 


= c 


s (d) 


= d 


o\ 


= a\a2 



b2bid2aibib2cdcd 



Now we fold the edges 02 and b 2 . 

p (a 1 ) = ax 

p (a 2 ) = b 2 

Po(6i) = &i 

Po(&a) = b 2 

Po(c) = c 

p (d) = d 



Finally, we compute the map #i : G 2 — > G. 



^i(oi) = a 

<7i(&i) = &cJc6 

9i(b 2 ) = cdcb 

gi(c) = bcdd 

gi(d) = ddcbd 

a 2 = a\b\b 2 d\b\b2cdcd 



2.2 Step 2: Triangulating annuli 

We can interpret the loops <jj as immersions a^: S 1 — > Gj. The preimage 
of the vertex set of Gj subdivides S* 1 into intervals, and the restriction of Uj 
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fold 




subdivision 




4 



final gluing 





fold 
subdivision 




Figure 3: Decomposition into annuli. 



to such an interval is a homeomorphism onto the interior of an edge in G{. 
Hence, we can label each interval with the corresponding edge in Gi. We 
refer to this construction as spelling er, along S l . 

Now, for each i G {0, . . . , 2n}, we take an annulus Ai and spell the word <7j 
along one boundary component and o"j + i along the other. We orient the two 
boundary components of Ai such that they are freely homotopic as oriented 
loops. 

This labeling defines a gluing of A, and Aj+i, and the homeomorphism 
9n '■ G 2n — G = G induces a gluing of A 2n -i and A Q (which we refer to as 
the final gluing), giving us the desired torus T (see Figure 0). Figure f| shows 
the gluing of A and A\ for Example |2.2| . 

Fix some i G {0, . . . , 2n} and spell ai along S 1 . Notice that each label on 
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c 



Figure 4: The gluing of Aq and A\ for Example [2.2| . 

S* 1 occurs twice, once for each direction (see Example |2.2| and Figure [|). If 
we identify corresponding intervals, we obtain the graph Gi, and if we take 
the cone over S 1 , remove the cone point, and identify intervals with identical 
labels, then we obtain the surface S ||BriOO| , Section 5.1]. 



Hence, we only need to extend the edge pairing on the boundary of the 
annulus Ai to an appropriate face pairing of a triangulation of all of A+. 
Recall that we want to choose the face pairing in such a way that the cone 
over T (with the cone point removed) becomes a 3-manifold when we identify 
corresponding triangles. 

We first find a suitable triangulation of an annulus A 2 i corresponding to 
a subdivision Sj. We decompose A 2 i into rectangles corresponding to edges 
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that are not subdivided and pentagons corresponding to those edges that are 
subdivided (see Figure £|). The edge pairing on the boundary of A 2i induces 
a pairing of rectangles (resp. pentagons). Finding a triangulation of A 2 i that 
is compatible with this pairing is straightforward. 

We now construct a suitable triangulation of an annulus A 2 i+i corre- 
sponding to a fold Pi. Edges that are not involved in the fold give rise to 
paired rectangles contained in A 2 i + i (see Figure £|), and as above, we easily 
find a triangulation of these rectangles that is compatible with the pairing. 

Let a, b denote the two edges involved in the fold, i.e., we have Pi(a) = 
Pi(b) = b'. By exchanging a and b or reversing the orientation of a and b 
as necessary, we may assume that the loop <7 2i+ i has a subpath of the form 
w = abua or w = abub, where u is a path that contains neither a nor b. For 
concreteness, we focus on the case w = abua. The first fold of Example |2]^ 
falls into this case, with w = a 2 b 2 b\a 2 . The construction in the remaining 
case is similar. 

The loop cr 2i+2 has a corresponding subpath of the form w' = u'b'. Let A 
be the triangle spanned by the initial endpoint of w' and the occurrence of a 
in w. We pair A with the triangle A' spanned by the terminal endpoint of 
w' and the occurrence of a in w (see Figure £|). 

Let Ai be the triangle spanned by the occurrence of b in w and the initial 
endpoint of w', and let A 2 be the triangle spanned by the occurrence of b' 
in w' and the terminal endpoint of u. Observe that after identifying paired 
edges, Ai and A 2 have a side in common, so we can think of Ai and A 2 as 
spanning a rectangle between b and b', which induces a triangulation of the 
rectangle spanned by the occurrence of b in <J 2 i+\ and the occurrence of b in 
o"2i+2 (see Figure £|). This completes the triangulation and face pairing of 
A 2 i + i, which completes our construction. 

Given the triangulation of A 2i+ i constructed above, we can think of the 
annulus A 2i+ i as interpolating between the graph G 2 i + i and the graph G 2 i +2 . 
In other words, we think of the fold as occurring continuously, by identifying 
larger and larger segments of the edges a and b (see Figure [|). Moreover, this 
continuous folding process is compatible with the embedding of the graphs 
in the surface S. This observation shows that the complex K and its trian- 
gulation have the desired properties, in particular Property [|. 

Summing up, we have obtained the following main result of this paper. 

Theorem 2.3. Let f be a homotopy equivalence of a finite graph G, repre- 
senting a homeomorphism <p of a once-punctured surface. Then the following 



10 



h b 2 



\ 

h b 2 



h b 2 



Figure 5: Continuous fold for Example |2.2| . 

is an effective procedure for computing a triangulation of the mapping torus 
ofo. 

1. Decompose the homotopy equivalence f : G —> G into a sequence of 
subdivisions and folds, followed by a homeomorphism (Section \2. 1\ ). 

2. Obtain the torus K as a gluing of one annulus for each subdivision and 
fold in the above decomposition (Figure^). 

3. Triangulate the individual annuli and construct a face pairing (Sec- 
tion \2~J\) . 

4- Construct a triangulation of the mapping torus of the surface homeo- 
morphism <p by taking the cone of K and glue tetrahedra according to 
the face pairing. 

□ 

The program jmt is an implementation of this procedure. 
2.3 Complexity analysis 

The purpose of this subsection is to obtain an estimate on the number of 
tetrahedra in the triangulation that we have constructed. 

Recall that we defined the size of a map h: G' — ► G to be the sum of 
the lengths (in the usual path metric) of the images of the edges of G' . For 
example, the size of the first map in Example ^| is 23. Let S(h) denote the 
size of h. 

We say that a map g : G' — ► G is tight if 

1. for every edge e of G', the restriction of g to the interior of e is an 
immersion, and 
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2. for every vertex v of G', there are two edges emanating from v that 
cannot be folded (not even after a subdivision). 



Note that tightness can always be achieved by homotopy. 

Proposition 2.4. Let f : G — > G be a tight homotopy equivalence repre- 
senting a homeomorphism <ft of a once-punctured surface of genus g, and 
assume that G has no vertices of valence less than three. Then the number 
of tetrahedra in the triangulation of is bounded by 16(5g — 2)S(f). 

Proof. Since folding reduces size and annuli come in pairs (a subdivision 
annulus followed by a folding annulus), the number of annuli is bounded 
above by 2S(f). A simple application of Euler characteristics shows that 
G has no more than 6g — 3 edges because the valence of each vertex is at 
least three. Subdividing and folding, however, may create additional edges 
as well as vertices of valence one or two, so we need to understand the effect 
of subdividing and folding on the number of edges. 

To this end, we introduce the notion of partial folds, i.e., folds where both 
participating edges have to be subdivided, and full folds, i.e., folds where at 
least one of the participating edges is not subdivided |BH92||. Clearly, a 



subdivision followed by a full fold does not increase the number of edges, 
whereas a subdivision followed by a partial fold increases the number of 
edges by one. 

A partial fold reduces the number of possible folds by one because the 
map resulting from it is an immersion around the new vertex created by the 
fold, so the only folds that are possible after a partial fold are those that 
were available before. 

Similarly, a full fold does not increase the number of possible folds. This 
means that the number of partial folds that occurs in our construction is 
bounded by the number of folds that the map / : G — > G admits. Since / is 
tight, the number of folds at one vertex v is bounded by val(f) —2. Summing 
up, we see that the number of possible folds is bounded by 

(val(^) - 2) = -2 X (G) = 4g-2. 

Hence, the number of edges after a fold is bounded by 6g — 3 + 4g — 2 = 
lOg — 5. Subdivisions increase the number of edges, so the number of edges 
at any point in our construction is bounded by lOg — 4. 
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The number of tetrahedra belonging to one annulus is bounded above 
by 8(5g — 2) (four tetrahedra per edge), which gives us a theoretical upper 
bound of 16(5(7 — 2)i>(/) on the number of tetrahedra in our triangulation of 
M . □ 

Similar arguments show that for fixed genus, the time it takes to compute 
a triangulation is linear in the size of the input. We note that in practice, 
partial folds seldom occur, and triangulations tend to be considerably smaller 



than the bound given in Proposition WA 



3 Solving the conjugacy problem 

The algorithm from Section |2| and SnapPea's isometry checker provide a 
practical way of testing a necessary condition for two pseudo-Anosov homeo- 
morphisms 0, 7 to be conjugate in the mapping class group. Namely, we can 
compute the mapping tori of and 7 as in Section §, and then SnapPea's 
isometry checker will determine whether the two mapping tori are isometric. 
If they are not isometric, we can immediately conclude that the and 7 are 
not conjugate. However, if the mapping tori are isometric, we cannot yet 
conclude that they are conjugate. The purpose of this section is to provide 
an effective sufficient criterion for conjugacy. 

Problem 3.1. Let 7 and be automorphisms of a surface S with one punc- 
ture. Both automorphisms are assumed to be presented as products of Dehn 
twists on S. 

The conjugacy problem asks for a decision procedure to determine whether 
or not 7 and are conjugate in the mapping class group of S. That is, does 
there exist an automorphism, ip, of S such that 7 = <y9 _1 0<^? 

The restricted conjugacy problem is slightly easier in that it assumes that 
is in fact pseudo-Anosov. As this is not a large restriction from now on we 
will assume that is pseudo-Anosov. 

A complete (albeit impractical) solution for the conjugacy problem has 
been given by Hemion |[Hem79|| . The restricted case has also been solved by 



Mosher ]Mo|8g . 
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3.1 Notation 



In order to describe our solution of the restricted conjugacy problem, we need 
to introduce some notation. Let 7 : S — > S and : S — > 5 be pseudo-Anosov 
homeomorphisms with isometric mapping tori. SnapPea will detect this and 
compute an isometry h: M 7 — > M^. 

Let i y : S — > M 7 and ^ be the two inclusion maps that realize 5 as a fiber 
of the induced fiber structures JF 7 and JF^, on M 7 and M^. Set F 7 = iy(S') 
and define in a similar fashion. 

Let p y : M 7 — > S 1 be the map from M 7 to the circle induced by JF 7 and 
define similarly. 

Let a G Mod(M0) denote a typical element of the isometry group of M^. 
As is hyperbolic, Mod(M^) is a finite group and can be computed by 
SnapPea. Set G a = a7i(F 7 ). 

Finally, pick any g G 7Ti(M 7 ) with the property that (p 7 )*<? = 1 G Z. We 
say that such a loop represents the S 1 - orientation. If (p^ahj^g equals 1 (—1) 
then we say that ah preserves (reverses) S 1 - orientation. 

3.2 Retriangulation and the fundamental group 

Before solving the restricted conjugacy problem we will need a pair of sub- 
routines which determine the images of elements of 7Ti(M 7 ) under the map 

(p<t,o-h)*. 

First, we need an algorithm which decides whether (p^,)*: ni(G a ) — > Z 
has trivial image. The idea here is to keep track of a set of generators for tti(S) 
under the maps i 7 , /i, cr and p^. Unfortunately, these homeomorphisms do not 
all respect common simplicial structures on M 7 and M^. To fix this problem 
one must find an appropriate set of generators in each retriangulation of 
M 7 . Then, once SnapPea finds a geometric triangulation of M 7 , we can push 
the generators onto the one-skeleton of the Ford domain. The isometry ah 
takes them to edge paths in the one-skeleton of where we can reverse the 
process. Finally, p^ projects the generators of 7r 1 (G , cr ) to the circle where it 
is easy to check whether or not they are all contractible. 

Second, we will need an algorithm which decides whether ah preserves 
S* 1 - orient at ion. To do this, construct any loop g G 7Ti(M 7 ) which represents 
the ^-orientation. As in the previous algorithm take the image of g under 
the map (p^ah)*. Check that this image is the positive generator of 7r 1 (S' 1 ). 

The bookkeeping problem of keeping track of surface subgroups of 7ii(M 3 ) 
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under retriangulation does not yet have an implemented solution. We would 
be very interested in the work of any reader who is willing to write such a 
program. It should be remarked that the subroutines above do a little more 
work than is strictly necessary. It would suffice to keep track of a two-chain 
representing the fiber of M 7 . Again, this is a straightforward problem which 
does not yet have an implemented solution. 



3.3 The algorithm 

Begin with two homeomorphisms 7,0: S — > S. If one of them fails to 



be pseudo-Anosov, then the software described in ||BriOO|| will detect this. 
Otherwise, construct the mapping tori of 7 and 0, M 7 and M$. If Snap- 
Pea reports that M 7 and are not isometric, we conclude that 7 cannot be 
conjugate to 0. This resolves the issue for a vast majority of possible pairs 
of 7 and 0. 

Now, suppose SnapPea reports the two mapping tori are isometric. We 
cannot yet conclude that the two automorphisms are conjugate. It may be 
that 7 and have homeomorphic mapping tori but are not conjugate because 
they give rise to distinct fiber structures in the resulting three manifold. Also, 
it may happen that 7 and induce identical fiber structures while reversing 
S 1 -orient at ion. In this case we show that 7 is conjugate to _1 . 

If 7r 1 (G CT ) has nontrivial image under (p^)* then clearly G a is not isotopic 
to F^. We conclude that if 7r 1 (G cr ) has nontrivial image in 717 (S 11 ) = Z for 
every a in Mod(M^) then 7 is not conjugate to 0. 

We claim that if there exists some a such that ah preserves ^-orientation 
and iTi(G a ) is contained in the kernel of the natural projection, then 7 is 
conjugate to 0, which completes our solution of the restricted conjugacy 
problem. The rest of this section is devoted to a proof of this claim. 



3.4 Isotoping G a 

Assume now that 77 (G>) has trivial image in (p^)* for some fixed o G 
Mod(M^). At this point we need a weak form of Theorem 4 from [ [Thu86fl : 



Theorem 3.2. G a is isotopic to a properly embedded surface that is either a 
leaf of or has only saddle singularities for the induced singular foliation 
ofG a . The boundary component of the isotoped G a is either a leaf of T^dM^ 
or is transverse to F^dM^. □ 
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We use this as follows: 



Corollary 3.3. Ifni{G a ) is in the kernel of (p<f,)* then G a is isotopic to Fa,. 

Proof. Suppose, to obtain a contradiction, that G a is not isotopic to Fa,. By 
Thurston's theorem we may isotope G a so that the induced foliation has only 
saddle singularities. 

Let Mi = Fa, x K. be the infinite cyclic cover of M coming from Ta,. By 
assumption we may lift G a to M%. Note that projection onto the second 
factor Mz — > K. gives a Morse function when restricted to G a . However this 
Morse function must have a maximum. The induced foliation of G a has a 
singularity at this maximum, but this singularity cannot possibly be a saddle 
singularity. This is a contradiction. 

□ 

Thus we may isotope ah so as to obtain G a = Fa,. Cutting along F 7 and 
Fa, we obtain a map h' : S x I —>■ S x I that takes S x and S x 1 to S x 
and 5x1, but not necessarily in that order. It may be that ah reverses the 
5 1 -orientation. 

It follows that h', and hence ah, is isotopic to a map that preserves fibers. 
(See, for example, Lemma 3.5 of |[Wal68|| .) Letting ho = (ah)\F^ we find that 



either 7 = h^ 1 4>h or 7 = h^ 1 0~ 1 /i o depending on whether ah preserves or 
reverses S^-orientation. This completes the proof of the claim and thus shows 
the correctness of our algorithm, which we sum up in the following theorem. 

Theorem 3.4. Let 7 and cf) be pseudo-Anosov homeomorphisms of a once 
punctured surface, S. The following is a procedure to decide whether the 
two mappings are conjugate in the mapping class group of S, if SnapPea is 
allowed as a subroutine. 



1. Apply Theorem 2.5 to obtain the mapping tori, M 7 and Ma, 



2. Using jsnap and SnapPea, determine whether or not M 7 and Ma, are 
isometric. If not, then 7 is not conjugate to </>. 

3. Using SnapPea, enumerate all isometries between M 1 and Ma,. 

4- Determine whether any of these are S 1 - orientation and fibre preserving. 
If none are, then 7 is not conjugate to 0. Otherwise, 7 is conjugate to 
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□ 



The complexity of this algorithm depends on the complexity of SnapPea, 
which we treat black box in this paper. 



A Generating input for SnapPea 

In order to generate input for SnapPea, the output of jmt has to be translated 
into SnapPea's triangulation file format by a second program (called jsnap). 
The purpose of this section is to discuss the intermediate format, which allows 
for quick and easy generation of input for SnapPea. We illustrate this format 
with a familiar example from |[lhu78[ | . 



The intermediate format admits two types of input lines, for tetrahe- 
dra and for gluings. An input line defining a tetrahedron has the form 
"T v\ v 2 t>3 t>4," where v±, .. . , v 4 are distinct labels of the vertices. Tetrahe- 
dra are glued implicitly if they have three vertices in common, and they can 
be glued explicitly by entering a line of the form "G v\ v% V3 wi W2 u>3," 
where t>i,t>2,t>3 and wi,W2,W3 are the labels of two faces of tetrahedra. In 
this gluing, the side [t>i,t>2] is glued to the side [wi,^], the side [^2,^3] is 
glued to the side [u^iite], etc. Empty lines and comments beginning with 
'//' are also allowed. 

Although the four vertex labels of a tetrahedron have to be distinct, two or 
more vertices of a tetrahedron may be identified after gluing. Distinct vertex 
labels are only needed in order to uniquely specify the sides of tetrahedra. 

Example A.l (Figure-eight knot). Consider the familiar figure-eight knot 
(see Figure An ideal triangulation of the complement of this knot can be 
expressed as a gluing of two tetrahedra [|Thu78| , Chapter 1] (see Figure 0). 
We encode the gluing as follows. 

T a b c d 
T b c d e 
G b e d a c d 
G c b e a b d 
G c e d a c b 

Feeding the above five lines into jsnap yields an encoding of the trian- 
gulation in SnapPea's file format. According to SnapPea, the fundamental 
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Figure 7: The figure-eight knot complement, expressed as a gluing of two 
ideal tetrahedra. Triangles are glued such that the line style and the direction 
of arrows are matched [[Thu78| , Chapter 1]. 



group of the resulting manifold M has the following presentation, which we 
modify by a sequence of Tietze transformations | LS77| , Chapter II. 2]. 

7Ti(M) = < x, y | yxxxyxyyx = 1 > (1) 

= < x, y , z | yxxxyxyyx = 1, y = zx > 

= < x, z I xzxxzxzxz = 1 > 

= < x, z I xzxzxxzxz = 1 > (2) 

Presentation (|2|) is the presentation of the fundamental group of the comple- 
ment of the figure-eight knot given in PZ85| , Example 3.8]. 



18 



Figure 8: A set of generators of the mapping class group. 



While SnapPea's format is an extremely efficient representation of trian- 
gulations, understanding it requires some effort. The program jsnap acts as 
an interface between SnapPea and the human user. If you can draw or visu- 
alize a triangulation of a 3-manifold, then you can also enter it into j snap. 



B Sample Computations 

We present some sample computations that illustrate the power of the soft- 
ware discussed here. We define surface automorphisms as compositions of 
Dehn twists with respect to the set of curves shown in Figure || We equip 
the surface with a normal vector field. When twisting with respect to a curve 
c, we turn right whenever we hit c (resp. left for inverse twists). The set of 
these Dehn twists generates the mapping class group |[Lic64| . 

We first present an example for which we can easily verify correctness. 

Example B.l (Figure-eight knot revisited). Let S be a punctured torus, 
and let 0: S — » S be given by 

<P = D a«D d0 . 

The complement of the figure-eight knot (see Figure |6|) is homeomorphic 
to the mapping torus of <fi (we will verify this soon). Given this composi- 



tion of Dehn twists, the train track software described in |PriOO|| determines 
that <p is pseudo-Anosov with growth rate A ~ 2.61803399, and SnapPea de- 
termines that the mapping torus M<f, is a hyperbolic 3-manifold of volume 
V « 2.02988321 with one torus cusp. 
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Moreover, the train track software computes the following topological 
representative / : G — > G of 4>, where G is a graph with one vertex and two 
edges. 

f(a) = ba 
fib) = bba 
a = abab 



This yields the following presentation of the fundamental group of the map- 
ping torus of (p. 

■K\{M^) = < a, b, t | tat = ba, ibt = bba > 

A few Tietze transformations show that 7Ti(M^) is the fundamental group of 
the figure-eight knot complement. 

tti(M^) = < a, b, t | tat = ba, ibt = bba > 
= < a, t | taataiaia = 1 > 
= < a, c | acacaacac = 1 > 



As in Example |A.1| , we have obtained the presentation given in [ |BZ85| , Ex 



ample 3.8]. Finally, the presentation of ttx(M^) computed by SnapPea is 

tti(M^) = < x, y | xyyyxyxxy = 1 >, 



which agrees with Presentation (|T|) in Example ATI 



Alternatively, we can run SnapPea's isometry checker on Example |A.l 



and Example [B.l| in order to see that we get the same hyperbolic 3-manifold 
in both examples. 



Example B.2 (Genus 3). Let S be a surface of genus 3 with one puncture, 
and let 0: S — > S be given by 

6 = D d D D, D r D H D- 1 . 

t do co ai ci 0,2 a2 

The train track software identifies as a pseudo-Anosov homeomorphism 
with growth rate A ~ 2.04249053, and SnapPea determines that the mapping 
torus is a hyperbolic 3-manifold of volume V ~ 4.93524268 with one torus 
cusp. 
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These applications only show a small part of all the possibilities. The 
train track software computes a plethora of information about surface home- 
omorphisms |BriOO|| , and SnapPea allows for a detailed analysis of (hyper- 



bolic) 3-manifolds. We believe that the combination of the two packages 
may become a valuable tool for topologists. 
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